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Abstract 

The model of strongly correlated electrons with the correlated hopping term 
and an additional interaction between holes V is solved exactly in one dimen- 
sion at a special point where the number of hole pairs is conserved. As a func- 
tion of the interaction V the phase diagram has a rich structure. There exist 
gapless phases with dominating density-density correlations (metallic) or pair- 
pair correlations (superconducting) as well insulating and phase-separated 
phases. The stiffness constant of the conductivity, the effective transport 
mass, and the compressibility are also exactly calculated in several regions of 
the phase diagram. 
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The study of strongly correlated fermion systems has attracted considerable attention 
during the last decade, in particular due to the discovery of high-temperature supercon- 
ductivity. Exactly solvable models, mainly those obtained using the Bethe-ansatz, have 
brought a very important progress in this area. Two well-studied models, which are exactly 
integrable in one dimension, are the Hubbard model and t — J model [1-5]. In [6] a new 
model was proposed, which is again solvable in one dimension, and it was shown [7] to be 
relevant to the theory of "hole" superconductivity introduced by Hirsh [8]. According to 
this theory the introduction of a correlated hopping term in the Hubbard Hamiltonian leads 
to an attractive effective interaction between the holes and the formation of Cooper pairs. 
Although this picture was confirmed by a BCS-type mean-field theory it is desirable to find 
exact results confirming this behavior. A first step in this direction has been made in [7] 
where a simplified version of Hirsch's Hamiltonian has been studied in one dimension. It 
was found that this model has a strong tendency towards superconductivity although it does 
not have finite off-diagonal long-range order. 

Another interesting general problem in the physics of highly correlated electronic systems 
is the characterization of the metallic, insulating and superconducting phases [9,10]. In this 
sense it is desirable to introduce exactly integrable models showing these phases. In this 
letter we introduce a model where all these phases are present. This model along with a 
correlated hopping term ti also has a static interaction V between holes. We show that this 
model is integrable provided that t\ = t, where t is the standard hopping matrix element, 
and study the superconducting properties at this point analytically. The phase diagram is 
very rich with a crossover from a region exhibiting long-range order, which corresponds to 
phase separation, to a regime with a complicated antiferromagnetic order, that corresponds 
to an insulating phase. Gapless regions with dominating superconducting pair correlation 
functions or dominating density- density correlations without long-range order also exist. 

The Hamiltonian of the correlated hopping model with an additional interaction beetween 
holes on a one-dimensional lattice with L sites and periodic boundary conditions is given by 
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H(t!,V)= -J2J2 [{ C ta C 3+l,<r + 4 +lt<T C jtff ) 
j=l a =±l 

X (l - tl^ + (l+ CT )/2 -a) 

+ V(2- n jja - n j+1>ff ) (1 - _„)], (1) 

where c Ji0 - is the annihilation operator for an electron with spin a at site j and n^ a = c^ c c^ a . 
Here we choose the hopping matrix element t — 1. The integrable correlated hopping chain 
[7] is a special case of (1) with V = 0. 

Firstly we will be interested in the special case H(V) = H(ti = 1,V). Note that the 
analogous construction have been considered for the Hubbard model in [10-12]. The results 
of this consideration are quite different from ours. 

In the Hamiltonian H(V) only one type of processes is allowed : a hole at the site (j, a) 
may hop to an occupied neighbour site (j ± 1, a) provided that the site (j + (a ± l)/2, — a) 
is empty. This process corresponds to the motion of the hole pair on the half of the lattice 
distance. In the Hamiltonian H(V) there is also the static interaction between holes which 
corresponds to an extremely anisotropic spin-spin interaction [12]. We consider hereafter 
the situation where we have only pair of holes [13]. In this case we consider hole pairs as 
"particles" and treat this problem within the framework of the Bethe-ansatz method. 

The discrete Bethe-ansatz equations are derived following the standard procedure by 
imposing periodic boundary conditions. Each state of the Hamiltonian is specifed by the 
set of charge rapidities pj representing the momenta of the hole pairs. All rapidities within 
a given set have to be different corresponding to Fermi statistics. These rapidities are 
determined by the Bethe-ansatz equations 

m 

e 2iL Pj = ( _ 1)m -l JJ exp [i{p . _ m) _ ie{p . jpi)] . 
1=1 

exp [-iQ(pj,pi)} = S(pj,pi)/S(pi,pj), 

Sfapi) = (l - + e ip '- +iw ) . (2) 

The energy and momentum of the corresponding state are given by 
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E = -2j2(cos Pj + V-fi)-2Lfi, P = (3) 

3=1 3=1 

where the chemical potential \i has been added to control the particle number. 

The above equations are valid for density of holes p = m/L < 2/3. For p > 2/3 we can 
choose as the reference state the state with unoccupied sites and treat again this problem by 
the Bethe-ansatz method. Each state of the Hamiltonian is now specified by the set of charge 
rapidities kj representing the momenta of the electrons. These rapidities are determined by 
the following Bethe-ansatz equations 

e"* = - ft exp£(fc, - kj) - iQ(kj, fc,)]. (4) 
i=i z 

The corresponding energy and momentum are given by 

n n 

E = -4LV + (4F - n - 2 cos kj) , P = Y, k i- ( 5 ) 

3=1 3=1 

We may consider (2,3) and (4,5) in the straighforward method [14]. In the thermodynamic 
limit (2) and (4) are replaced by the integral equations for the distribution functions 

2nR(U) = K^U) (l - - J U ° K 2 {U - U')R(U')dU', 

Uo [\p, P< 2/3, 

R(U)dU = { (6) 
~ Uo I(l-P), 2/3<p<l, 

where for — 1 < V = — cos 7 < 1 

K a {U) = sin(o;7)[cosh U — cos(a7)] _1 , (a = 1, 2) 
and for V = — cosh A < — 1 

K a (U) = sinh(aA)[cosh U — cos(aA)] -1 , (a = 1, 2). 
The energy is given by 

2l B = C e{U)dU (?) 
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where the dressed energy e(U) is the solution of the integral equation 

1 r u ° 

e(U) = e (U) - — / K 2 (U - U')e(U')dU'. (8) 

271 J-Uo 

Here for V < 1 

eo(U) = -2sin 7 K 1 (f/) + ^(1 ± 3)(/i + 2 cos 7 ), 
with sign + (-) for p < 2/3 (p > 2/3), respectively. For V < -1 

e (C/) = -2 sinh AXi(C/) +2(^ + 2 cosh A). 

In all cases e(±C/ ) = 0. 

At 1/ < — 1 and p > 2/3, there is a gap in the energy spectrum and it corresponds to the 
existence of long-range order in the system. In this case the equations (6,8) are not valid. 

At V > 1 there is a gap for the arbitrary concentration of holes and this state has 
long-range order which corresponds to phase separation at V = — 1. 

For — 1 < V < 1 there is no gap in the eigenspectrum and in order to understand the 
superconducting properties of the model we shall investigate the long-distance behavior of 
the correlation functions. For this purpose we shall use two-dimensional conformal field 
theory [15] and analytic methods [16,17] to extract finite-size corrections from the Bethe- 
ansatz equations. The results of these calculations are the following. 

The long-distance power-law behavior of the density- density correlation functions is given 

by 

(p(r)p(O)) ~ p 2 + A x r~ 2 + A 2 r~ a cos{2k F r), (9) 

where 

2k F = Tip, p(r) = J2 CraCra- 

a 

The critical exponent a is expressed in terms of the dressed charge £(Uo) which satisfies the 
integral equation 
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1 r u o 1 
£(U) + — / K 2 (U - U')t{U')dU> =l--p, (10) 

The Eq. (10) and (14) add to the integral equations (6) and (8) [18]. The superconducting 
properties of the system manifest themselves in the behavior of the pair correlation function 

G P {r) = (c+ T c+ iCoi c OT ) ~ Br~ p , 

The exponent (3 is plotted in Fig.l for some values of the parameter V. These results 
are obtained by solving numerically the integral equations (6) and (10). We see from this 
figure that for all values of V there exists a critical density p c separating a regime with 
dominant density-density correlations (5 > 1 from a regime with dominant pair supercon- 
ducting correlations (3 < 1. These two regimes are represented in the phase diagram of Fig. 
2 by the phases C and B, respectively. The critical curve p c (V) is obtained numerically 
and is represented by the broken line. In phase B, where the superconducting correlation 
function decays slower than the density-density correlation function we have the closest 
analogy with a true superconducting phase. Analogous behavior of correlation functions 
are observed in other models [7,19], but differently from these models, in the present case 
we have the nice feature that for increasing values of p, the exponent (5 tends toward zero. 
This means that the system try to have a true off-diagonal long-range-superconducting order 
(ODLRSO). However at the same time when (3 — > 0, which happens for V — > 1~, the hole 
compressibility, as we shall see, tends to infinity and phase separation takes place. Then at 
V = 1 we have two phase transitions happening simultaneously. We may expect that in two 
dimensions these phase transitions happen at different places [20]. 

In region D (V < —l,p > 2/3) we have no bound pairs but a long-range "antiferromag- 
netic order", where hole pairs and electrons are located alternately. This is an insulating 
gapped phase. In order to see this we calculated the electrical conductivity in our model. 
This can be done by studying the Hamiltonian (1) with twisted boundary condition, where 
the twisting angle corresponds to an enclosed magnetic flux in the ring [9]. The conductiv- 
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ity is directly proportional to the charge stifness D c , which can be obtained from the change 
AE = D c (f) 2 /L of the ground-state energy E , for small and large L. From the charge 
stiffness we can define an effective transport mass by the relation 

m/m e = D°JD C , (12) 

where D® = 2sin(7rp)/7r is the charge stiffness of the non-interacting system and m e is the 
electron mass. Extending the Bethe ansatz equations (2-5) for the case of twisted boundary 
conditions [21] cf +lcr = exp(±i0 (7 )c^ (T with + = 0_ = and exploring the finite size 
corrections obtained from conformal invariance [15] we obtain, after some calculation 

D c = 2v F e/n, v F = e'(U )/(A7rR(U )), (13) 

where by prime we mean derivative. In Fig. 3 we show m/m e obtained numerically for 
several values of V . We clearly see in this figure that for V < —1 as p — > 2/3 the transport 
mass diverges, which indicates that phase D in Fig. 2 is an insulating phase, and at p = 2/3 
we have a Mott transition. The anomalous behavior in Fig. 1 and 3 around p = 2/3 
at V = —0.98 is due to the strong competition between the repulsive and kinetic energy, 
which for larger values of V will lead to the Mott transition. This figure also shows that for 
— 1 < V < 1 the carrier mass is around 2m e for moderate densities, which means that the 
conductivity is ruled by the motion of pairs. 

For V > 1 (phase A in Fig. 2) we have bound pairs, for any density and there is 
a long-range order in the system. This is a region where phase separation takes place. 
The hole pairs are separated from the particles. In order to see this we calculated the hole 
compressibility k — 1/ (p 2 ^¥-), where E p is the ground-state energy for a given concentration 
p. Exploring, as before, the finite-size corrections obtained from conformal invariance [15] 
we obtain k = D c /2v F p 2 . In Fig. 4 we show the curves of k for certain values of V. We 
clearly see in this figure that k increases drastically for all densities as we tend toward V = 1, 
indicating a phase separation for V > 1 (phase A in Fig. 2) [22]. 

Thus we have investigated a correlated hopping model [7] with an additional hole-hole 
interaction V at the special point t± — 1 where the number of hole pairs is conserved and 
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we have only hole pairs. We have shown that our model exhibits phase separation and 
Mott transition at V = 1 and V = — 1, respectively. At — 1 < V < 1 there are regions 
with dominating superconducting correlations. Although the system in this region has no 
ODLRSO it is very close to the true superconducting state since the critical exponent of the 
superconducting correlation function (5 tends to zero. 

In the general case at < t± < 1 desintegration of hole pairs is possible, and the number 
of hole pairs is not conserved. However these states are not favorable in energy and their 
contribution is negligible since they are separated from those considered above by a large 
energy gap, at least for a small density of holes. As a result, we may expect that at arbitrary 
ti the qualitative behavior of the system will be same as that considered at t± — 1. Certainly 
it is desirable to investigate this problem using perturbation theory and exact diagonalization 
of small systems. 

This work was supported in part by Conselho Nacional de Desenvolvimento Cientifico 
e Tecnologico - CNPq - Brazil, and by Fundagao de Amparo a Pesquisa do Estado de Sao 
Paulo - FAPESP - Brazil. 
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FIGURES 

FIG. 1. The exponent (3 of the pair correlation function as a function of hole density p for some 
values of the interaction V. 

FIG. 2. Phase diagram as a function of the interaction v and hole density p. Phases A and 
D are gapped and correspond to phase separations and insulating phases, respectively. Phases 
B and C are gapless and corresponds to the region with dominating behavior of superconducting 
correlations and density-density correlations, respectively. 

FIG. 3. The effective transport mass as a function of density of hole p for some values of V. 

FIG. 4. The hole compressibility as a function of density of hole p for some values of V . 
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